Elastic Behavior of a Two-dimensional Crystal near Melting 
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Using positional data from video-microscopy we determine the elastic moduli of two-dimensional 
colloidal crystals as a function of temperature. The moduli are extracted from the wave-vector- 
dependent normal mode spring constants in the limit q — > and are compared to the renormalized 
Young's modulus of the KTHNY theory. An essential element of this theory is the universal predic- 
tion that Young's modulus must approach 167T at the melting temperature. This is indeed observed 
in our experiment. 
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In the early 70th Kosterlitz and Thouless Q devel- 
oped a theory of melting for two-dimensional systems. 
In their model the phase transition from a system with 
quasi-long-range order |2| is mediated by the unbinding 
of topological defects like vortices or dislocation pairs in 
the case of 2D crystals. They showed that the phase 
with higher symmetry has short-range translational or- 
der. Halperin and Nelson 0, 0] pointed out that this 
phase still exhibits quasi-long-range orientational order 
and proposed a second phase transition now mediated 
by the unbinding of disclinations to an isotropic liquid. 
The intermediate phase is called hexatic phase. This the- 
ory, being based also on the work of Young Q , is known 
as KTHNY theory (Kosterlitz, Thouless, Halperin, Nel- 
son and Young); it describes the temperature-dependent 
behavior of the elastic constants, the correlation lengths, 
the specific heat and the structure factor (for a review see 
0,0)- Experiments with electrons on helium B, and 
with 2D interfacial colloidal systems EH, El Il2l fl3 . 0] 
as well as computer simulations |l5l Il7| have been 
performed to test the essential elements of this theory. 
But research has mainly focused on the behavior of the 
correlation functions (an illustrative example is the work 
of Murray and van Winkle 11]). Only a few works can 
be found that deal with the elastic constants, especially 
the shear modulus, 0,0,0,113 even though the Lame 
coefficients and their renormalization near melting take 
a central place in the KTHNY theory. 

A very strong prediction of the KTHNY theory has 
never been verified experimentally. It states that the 
renormalized Young's Modulus K R (T), being related just 
to the renormalized Lame coefficients ~p, R and Xr, must 
approach the value 16-7T at the melting temperature 0, 

K R {T) = 47^(1 - Ji R /(2p R + X R )) -> 16tt if T — > T m , 

(1) 

which is obviously an universal property of 2D systems at 
the melting transition. This Letter presents experimen- 
tal data for elastic moduli of a two-dimensional colloidal 
model system, ranging from deep in the crystalline phase 



via the hexatic to the fluid phase. These data, indeed, 
confirm the theoretical prediction expressed by eq. ljT]l. 

The experimental setup is the same as already de- 
scribed in [23|. The system is known to be an al- 
most perfect 2D system; it has been successfully tested, 
and explored in great detail, in a number of studies 
lid l2ll \'2'A |23| . Therefore we only briefly summarize the 
essentials here: Spherical colloids (diameter d = 4.5 /jm) 
are confined by gravity to a water/air interface formed 
by a water drop suspended by surface tension in a top 
sealed cylindrical hole of a glass plate. The flatness of 
the interface can be controlled within ± half a micron. 
The field of view was 835 x 620 pm 2 containing typi- 
cally up to 3 • 10 3 particles (the whole system has a size 



of 50 



and contains about 3 • 10 particles). The 



particles are super-paramagnetic, so a magnetic field B 
applied perpendicular to the air/water interface induces 
in each particle a magnetic moment M = \B which 
leads to a repulsive dipole-dipole pair-interaction energy 
of (3v(r) — r/(y / 7r/5r) 3 with the dimensionless interaction 
strength given by T = f3(fi /4:Tr){xB) 2 (irp) 3 / 2 (J3 = 1/kT 
inverse temperature, x susceptibility, p area density). 
The interaction can be externally controlled by means 
of the magnetic field B. T was determined as in Ref. 0] 
and is the only parameter controlling the phase-behavior 
of the system. It may be considered as an inverse reduced 
temperature, T = 1/T. For V > T m = 60 the sample is a 
hexagonal crystal |14i |22j . Coordinates of all particles at 
equal time steps and for different 'temperatures', i.e. T's, 
were recorded using digital video-microscopy and evalu- 
ated with an image-processing software. We measured 
over 2-3 hours and recorded trajectories of about 2000 
particles in up to 3600 configurations, for a large number 
of different T's ranging between T = 49, deep in the fluid 
phase, to r = 175 in the solid phase. These trajectories 
were then further processed to compute the elastic con- 
stants of the colloidal crystal as a function of the inverse 
temperature T. 

Our data analysis is based on the classical paper of Nel- 
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son and Halperin (NH) on dislocation-mediated melting 
in 2D systems |4|. Their considerations start from the 
reduced elastic Hamiltonian, 



mi 



(2) 



where a is the lattice constant of a triangular lattice (next 
neighbor distance), while ~p = fia 2 /kT and A = Xa 2 /kT 
denote the dimensionless Lame coefficients (/3 = 1/kT). 
Uij(r) = (d rj Ui(f) + d ri Uj(f))/2 is the usual strain ten- 
sor related to the displacement field u(r). At tempera- 
tures T near the melting temperature T m the field Uij if) 
contains singular parts u^j n9 (r) due to dislocations; it 
can be decomposed into Uij (r) = u*j™ 9 (r) + 4>ij (?) with 
i'ijif) = (9 rj 4>i(f) + dr i 4>j(r))/2 being a smoothly vary- 
ing function (<j)(r) is the regular part of the displace- 
ment field u(r)). When this decomposition is inserted 
in eq. the Hamiltonian decomposes into two parts, 
I3H E = 1/2 J d 2 r/a 2 {2jl4>l +\4>l k )+/3H D with H D rep- 
resenting the extra elastic energy that is due to the dis- 
locations. NH were able to derive a set of differential 
equations for renormalized Lame coefficients, ~p R X R , by 
means of which He can again be written as in eq. (J2J, 



mi 



d 2 r 



[ZpRU>ij + X R u 2 kk ] 



(3) 



Because the effect of the dislocations are entirely ab- 
sorbed into the elastic constants, the strain tensor in 
eq. ® can now be assumed to be again regular every- 
where and for all T < T m . 

Our experiment measures the trajectories Ti{t) of N 
particles of a colloidal crystal over a finite time window of 
width t exp . Associating the average (ri)t cnsp with a lattice 
site Ri, we can, for each particle, compute displacement 
vectors u(Ri) = ri — Ri . The Fourier transforms of these 
displacement vectors, u(q) = N^ 1 / 2 ^^e l ^ R u(R), are 
now used for the numerical computation of renormalized 
elastic constants. This has been done in the following 
way. Starting from eq. (3.29) of the NH paper, 



kT 

lim iq 2 (u*(q) Uj (q)) = 

q^Q Vq 



—5i 



1-iR + X R qiqj 



fj, R (2^ R + X R ) q 2 



J4) 

we find, after decomposing the displacement field u(q) 
into parts u\ \ (q) and u± (q) , parallel and perpendicular 
to q, that 

^M+M = l^C,,^, C||(^:=( 9 2 (l«||( 9 1| 2 »lfe) 

! <K(<z)| 2 »(^ 



kT 



q^O 

iimC±(9), CAq) ■= {q' 



where vq = \/3a 2 /2 is the area per colloid in a triangular 
lattice. 



These two equations are the central equations for our 
data evaluation scheme, and are therefore more care- 
fully discussed. Deep in the solid phase (r)&t converges 
with increasing measurement time At to lattice sites, 
and (|m(.R)| 2 ) and (g 2 (|t?(<7)| 2 )) g _ t o remain finite, and lead 
thus to two non-zero elastic moduli: the shear modulus 
fi and the bulk modulus B = X + fi. By contrast, in the 
fluid phase (r)At will neither converge to, nor correlate 
with, any sort of lattice site and, in addition, as the mean 
square displacement is unbound, (<7 2 (|t?((7)| 2 )) g _ + o — * oo, 
so, both moduli will vanish (while B survives |2J|). 

Though reasonable at first glance, this interpretation 
of eqs. H5I6J) is in fact an oversimplification and ignores 
the limited applicability of both equations. This can best 
be seen by re-deriving them, first by Fourier-transforming 
eq. Q , by then inserting the decomposition u = u\ | + u± 
and by applying finally the equipartition theorem. The 
intimate relationship between eq. J3J and eqs. (TiHil) let 
us realize that the u(q) in eqs. (|5I6(I refer to a coarse- 
grained and thus regular displacement field, just as in 
eq. J2Jl. In other words, with eqs. il~>l(ill the softening of 
the elastic constants for T — > T,~ is inferred indirectly, 
namely from the change of behavior of the regular parts 
of a coarse-grained displacement field. 

We here identify this coarse-grained displacement field 
with the u(f(t)) = (r)t ~'> f (t) evaluated from our exper- 
imental data. In doing so, we have to be aware of the fol- 
lowing two points, (i) In an imperfect crystal, especially 
in the presence of dislocations, (r) tc;rp _ >00 does not always 
converge to lattice sites. As we are interested in the limit 
q — > 0, this is unproblematic, as long as these extra sites 
do not move. Looking at our measured trajectories, we 
have observed only thermally activated dislocation pairs, 
but no static, isolated dislocations traveling through the 
crystal, (ii) The displacement field computed from the 
experimental data has (and must have) parts stemming 
from dislocations. The resulting error should be small be- 
low T m when the number of dislocations is still relatively 
low (even at T = T m , the probability that a particle be- 
longs to a dislocation is only 1 % !). However, the error 
should become appreciable at T > T m where eqs. <|5I6[) 
can not be expected to strictly hold any longer. 

For the pair-potential ~ T/r 3 the elastic constants can 
be calculated in the limit r — > oo (T = 0) using sim- 
ple thermodynamical relations involving essentially lat- 
tice sums of the pair-potential. One finds X + ~p = 3.46 T 
and Ji = 0.346 T (2||. For convenience, we divide in the 
following all moduli by V. Fig. ((TJ shows (A + 2^)u /rfcT 
and iivo/TkT, obtained from this T = calculation, as 
thick solid arrows, and compares it to the expressions 
0i(q)/T and (±(q)/T from eqs. (|5I6() . as obtained from 
the measured trajectories for three different values of T. 
Let us first focus on the measurement for T = 75, 125 
where the system is deep enough in the crystalline phase 
for the assumption T = to be valid. C||(<7)/r and 
(±(9)/!^ indeed, tend to the predicted elastic constants 




FIG. 1: C||(<zV r in a ) and Cx(g)/T in b), as denned in 
cqs. (ITiK it . Each quantity is plotted for two different directions 
in q-space (F — > M and F — » if in the first Brillouin zone), 
and for three different values of the interaction-strength pa- 
rameter V: r — 52 in the fluid phase (open triangles), r = 75 
(open squares), and V = 125 (filled circles) in the crystalline 
phase. Thick solid arrows for a T = prediction of the elas- 
tic moduli, dashed solid lines for the predictions of harmonic 
crystal theory. For r = 52 just one band is shown. 



in the limit q — > 0, in agreement with eqs. ll">l(il) . For 
either of the two plotted quantities, we consider two dif- 
ferent high-symmetry directions in q space, which are 
<fi = qbi and q\ = q(b\ + ^Vv^ where b\ = (l,0)a and 
&2 = (1, \/3)a/2 are basis vectors of the (hexagonal) re- 
ciprocal lattice. At wavelengths larger than the lattice 
constant (qa < 1), the results for both bands are identi- 
cal, thus indicating an essentially isotropic u{q) at small 
q. C,\\{q) and C±(q) can also be associated with the q- 
dependent normal-mode spring constants (elastic disper- 
sion curves) of the discrete crystal, having longitudinal 
^longio) and transversal A tra „ s (q) branches. This can be 
(and have been 23]) compared to the band-structure pre- 
dicted by harmonic lattice theory (thick dashed lines in 
Fig. QJ, for details see H^|). In other words, what we do 
here is to derive elastic constants from the q — > behav- 
ior of the elastic dispersions curves (lim^o Xi ong (q) = 
(A + 2^)q 2 v Q , lim 9 ^ ^trans(q) = vq 2 v ). 

While for all our measurements above r m = 60 the 
resulting bands lie on top of the dashed thick lines in 
Fig. (JJJ , one finds a systematic shift to smaller values for 
T < r m . Fig. 0} shows, as an example, one out of the 
four bands (belonging to £j_(<f)/r) of the measurement 
in the fluid phase (r = 52). It lies an order of magnitude 
below the crystalline bands. 

In order to infer functions hr(T) and Afl(r) from these 
bands, we need to take the limit q — ► 0. Since at low q 
we have to expect finite size effects, and at high q, near 
the edges of the first Brillouin zone, effects resulting from 
the band dispersion of the discrete lattice, we choose an 
intermediate q regime (0.8 < qa < 2.5), indicated by the 
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FIG. 2: Elastic moduli of a 2D colloidal crystal as a function 
of the inverse temperature, obtained from extrapolating the 
bands in Fig. Q down to q = 0. The melting temperature is 
at r m = 60. Thick dashed lines for a T — prediction of the 
elastic moduli, A + 271 = 3.806r and 71 = 0.346F ;25j; thick 
solid line for the theoretical elastic constants, renormalized as 
described in Q. 

thick solid bar in Fig. 0), to extrapolate the bands to 
q = 0, applying a linear regression scheme. The extrap- 
olation procedure was optimized at high T for which we 
have a precise idea what constants we should find. For 
each modulus, extrapolation of the two bands depicted in 
Fig. 0} were checked: while for Cl|A\ one finds for both 
bands the same modulus, the upper band of the two for 
C_i_/T gave much better result and was henceforth taken. 
Also, different extrapolation schemes have been checked, 
but linear extrapolation turned out to produce a tolera- 
bly small error, much smaller than the main error of our 
measurement, estimated here from the standard error in 
the linear regression scheme. 

Fig. (0 shows the resulting moduli, for all values of 
r studied. Black symbols refer to systems in the crys- 
talline state (r > r m ), grey data points to those in 
the fluid/hexatic phase. We postpone the discussion of 
the data points at T < T m and first concentrate on the 
crystalline regime where a renormalization of Lame co- 
efficients really makes sense. The thick dashed lines in 
Fig. J2J represent the T = calculation which holds down 
to r values close to T — 75. The thick solid line shows the 
theoretical curve for /ir(T) and nn(r)+2Xjt(T), which we 
computed following the renormalization procedure out- 
lined in the NH paper (Eq. (2.42), (2.43) and (2.45) in 
with K(T) = 1.258 T at T = and T m = 60 as input 
to set up the boundary conditions.) Theory and exper- 
iment agree well, considering that no fit parameter has 
been used. For T > T m , all our results are converged, 
meaning that the computed moduli do not depend on the 
length of the trajectory. This is demonstrated by means 
of the r = 100 measurement for which the moduli were 
computed taking only the first third of all configurations 
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FIG. 3: Young's modulus, eq. (0, as a function of the inverse 
temperature, evaluated with the experimentally determined 
Lame coefficients of Fig. @ (symbols). The solid curve is 
Kr(T) renormalized according to Q, while the dashed curve 
is based on the T=0 prediction. 



(open square symbol in Fig. J2J). 

Fig. now checks eq. (Q, with Kr(T) evaluated us- 
ing the elastic moduli from Fig. J2J. Using the theoretical 
values from the T = calculation, we obtain K(T) = 
1.258 r, shown in Fig. @ as dashed line. The thick 
solid line shows the theoretical curve for Kr(T) which 
we computed with Lame coefficients that were renormal- 
ized following the NH procedure explained above. The 
main result of this work is that the experimental data 
points closely follow the theoretical curve and indeed, 
they cross 167r at T = T m in excellent agreement with 
the predictions of NH. The length of the remaining error 
bars correlate with the total measurement time, and the 
g-range chosen in the extrapolation step. 

The data points for T < T m in Fig. {Q and © should 
be treated with utmost caution. For the reasons given 
above, it is not clear to us whether or not eqs. Ij5l6|l is at 
all meaningful outside the crystalline phase. But even if 
it were, we should be aware that the results will always 
depend on the measurement time. This is demonstrated 
for a system at T = 49 for which the moduli were cal- 
culated taking again just a fraction of 1/3 of all config- 
urations (open symbols in Fig. J2J0- There is almost an 
order of magnitude difference to the data points based on 
all configurations, thus indicating the dependence of the 
moduli on the length of the analyzed trajectories. Phys- 
ically, one could interpret this in terms of a frequency- 
dependent shear modulus which for non-zero uj ~ {/t exp 
is known to exist even in fluids. 

To conclude, we have measured particle trajectories of 
a two-dimensional colloidal model system and computed 
elastic dispersion curves which at low q give access to the 
elastic constants. We thus measured /Ur, A# and Young's 



modulus Kr as a function of the inverse temperature 
r. All three quantities compare well with corresponding 
predictions of the KTNHY theory. Young's modulus, 
in particular, tends to 167T when the crystal melts, as 
predicted in |4J. 
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